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Using the relationship between the two-particle overlap functions (TOF's) and the two-body 
density matrix (TDM), the TOF's for the 16 0(e, e'pp) 14 C g . s . reaction are calculated on the basis 
of a TDM obtained within the Jastrow correlation method. The main contributions of the removal 
of 1 So and 3 Pi pp pairs from 16 O are considered in the calculation of the cross section of the 
ls O(e, e'pp) 14 Cg. s . reaction using the Jastrow TOF's which include short-range correlations (SRC). 
The results are compared with the cross sections calculated with different theoretical treatments of 
"pg"' the TOF's. 
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. Since a long time electromagnetically induced two-nucleon knockout has been devised as the most direct tool to 
study the properties of nucleon pairs within nuclei at short distance and thus the dynamical SRC in a nucleus |IJ ^) . 

Two nucleons can be naturally ejected by two-body currents, which effectively take into account the influence of 
subnuclear degrees of freedom like mesons and isobars. Direct insight into SRC can be obtained from the process 
where the real or virtual photon hits, through a one-body current, either nucleon of a correlated pair and both 
nucleons are then ejected from the nucleus. The role and relevance of these two competing processes can be different 
in different reactions and kinematics. It is thus possible to envisage situations where either process is dominant and 
. various specific effects can be disentangled and separately investigated. 

Various theoretical models for cross section calculations have been developed in recent years in order to explore the 
(S| ! effects of ground-state NN correlations on (e, e'NN) § |, |, |, @, §) and (7, NN) |[ [h], || || [H |1 @ knockout 
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PACS numbers: 21.60.-n, 21.10.Jx, 25.30.Dh, 27.20.+n, 24.10.-i 



I. INTRODUCTION 



1 reactions. It appears from these studies that the most promising tool for investigating SRC in nuclei is represented 



by the (e,e'pp) reaction, where the effect of the two-body currents is less dominant as compared to the {e,e'pn) and 
processes. Measurements of the exclusive 16 0(e, e'pp) 14 C reaction performed at NIKHEF in Amsterdam 
18| ] and MAMI in Mainz |T^, have confirmed, in comparison with the theoretical results, the validity of the 
direct knockout mechanism for transitions to low-lying states of the residual nucleus and have given clear evidence 
of SRC for the transition to the ground state of 14 C. This result opens up good perspectives that further theoretical 
and experimental efforts on two-nucleon knockout will be able to determine SRC. 

One of the main ingredients in the transition matrix elements of exclusive two-nucleon knockout reactions is the 
two-nucleon overlap function. The TOF contains information on nuclear structure and correlations and allows one 
to write the cross section in terms of the two-hole spectral function J2|. The TOF's and their properties are widely 
reviewed, e.g., in [ pl| . In an inclusive reaction, integrating the spectral function over the whole energy spectrum 
produces the TDM. 

In H the TOF's for the 16 0(e, e'pp) 14 C reaction are given by the product of a coupled and fully antisymmetrized 
pair function of the shell model and a Jastrow- type correlation function which incorporates SRC. Only the central 
term of the correlation function is retained in the calculation. 

A more sophisticated treatment is used in 0, where the TOF's are obtained from an explicit calculation of the 
two-proton spectral function of 16 p^ j, which includes, with some approximations but consistently, both SRC and 
long-range correlations (LRC). The numerical predictions of this model are in reasonable and in some cases in good 
agreement with data (l7|, [0| ML 

Although satisfactory, these first comparisons with data have also raised problems that require further theoretical 
investigation and a more refined treatment of nuclear structure aspects in the calculation of the spectral function. 
The study of the two-hole spectral functions including different types of correlations, however, requires substantial 
efforts in computational many-body physics and represents a very difficult task. 

A different method to calculate the TOF's has been suggested in [^3| using the established general relationships 
connecting TOF's with the ground state TDM. The procedure is based on the asymptotic properties of the TOF's 
in coordinate space, when the distance between two of the particles and the center of mass of the remaining core 
becomes very large. This procedure can be considered as an extension of the method suggested in j24|, where the 
relationship between the one-body density matrix and the one-nucleon overlap function is established. The latter 
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has been applied [|25|, |26|, 27, |28|, |29|, |30[ BlL |32| to calculate the one-nucleon overlap functions, spectroscopic factors 
and to make consistent calculations of the cross sections of one-nucleon removal reactions (p,d), (e,e'p), and (7,p) 
H ||, ||, H, H on 16 ||, |(| H and 40 Ca @ ||. Various correlation methods, such as the Jastrow method, 
the Green function method, the correlated basis function method and the generator coordinate method have been 
used to obtain the one-nucleon overlap functions which are necessary for cross section calculations. 

The first aim of the present paper is to apply the procedure suggested in [^3| to calculate TOF's for 16 O using the 
TDM calculated in |3j| with the Jastrow correlation method (JCM), which incorporates the nucleon-nucleon SRC. 
As a second aim, the resulting two-proton overlap functions are used to calculate the cross section of the 16 0(e, e'pp) 
reaction for the transition to the ground state of 14 C. The cross sections are calculated on the basis of the theoretical 
approach developed in || |[ |?J . The choice of the Jastrow TDM is determined by the convenience of its analytical 
form obtained in |Q, which makes practically possible the calculation of the TOF's. Of course, the reliability of the 
TOF's obtained in our method depends strongly on the availability of realistic TDM's. So, the usage in our work of 
the Jastrow TDM, though incorporating only SRC (and using harmonic-oscillator single-particle wave functions in 
the Slater determinant), must be considered as a first attempt to use an approach which fulfils the general necessity 
the TOF's to be extracted from the TDM and to apply them to cross section calculations of two-nucleon knockout 
reactions. 

The method to calculate the TOF's on the basis of the TDM is briefly outlined in Section II. The results of the 
calculations of the TOF's and the cross section of the 16 0(e, e'pp) 14 C g . s . reaction are presented and discussed in 
Section III. Some conclusions are drawn in Section IV. 

II. TWO-BODY DENSITY MATRIX AND OVERLAP FUNCTIONS 

In this Section we present shortly the definitions and some properties of the TDM and related quantities in both 
natural orbital (geminal) and overlap function representations. The method to extract the TOF's from the TDM 
used in this work is also given. 

The TDM is defined in coordinate space as: 

p^(x 1 ,x 2 ;x[,x / 2 ) = (^\a^x 1 )a^x 2 )a(x' 2 )a(x' 1 )\^), (1) 

where I'J'^) is the antisymmetric ^4-fermion ground state wave function normalized to unity and a^(x), a(x) are 
creation and annihilation operators at position x. The coordinate x includes the spatial coordinate r and spin and 
isospin variables. The TDM p^ is trace-normalized to the number of pairs of particles: 

TrpW =lj p^(x 1 ,x 2 )dx 1 dx 2 = A(A ~ 1] . (2) 

Since p' 2 ' is a Hermitian matrix its eigenstates ipa^ form a complete orthonormal set in terms of which p^ can be 
decomposed as 

pW(x 1 ,x 2 ;x' 11 x' 2 ) ^^Ai 2 )^ 2 )*^,^)^^!,^) • (3) 

a 

The eigenfunctions ip„ '(xi,x 2 ) are called natural geminals and the associated real eigenvalues Xa - natural geminal 

occupation numbers pi] . As a consequence of the antisymmetry of the nuclear ground state, the eigenvalues A^ 2 ' 
obey the inequalities: 



< AL 2) <(A- l)/2 for A odd, 
< Ai 2) < A/2 for A even. 



(4) 



The upper bound in Eq. (Q) is actually reached only for systems which are maximally correlated, as, e.g., the 
occupation number of zero-coupled pairs in the seniority formalism in the limit of large shell degeneracy. 

Of direct physical interest is the decomposition of the TDM in terms of the overlap functions between the A-particle 
ground state and the eigenstates of the (A — 2)-particle systems, since TOF's can be probed in exclusive knockout 
reactions. 

The TOF's are defined as the overlap between the ground state of the target nucleus vff^) and a specific state 
of the residual nucleus (C = A — 2) |2l| : 

$ a {x u x 2 ) = (*{p)|o(a;i)o(a: 2 )|*W). (5) 
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Inserting a complete set of (A — 2) eigenstates \a(A — 2)) into Eq. (|l|) one gets 

p {2) (x 1 ,x 2 ;x' 1 ,x' 2 ) =Y^K{xi,x 2 )^ a {x' 1 ,x' 2 ). (6) 

a 

The norm of the two-body overlap functions defines the spectroscopic factors 

S<?> = <**l*a>. (7) 

As in the case of the single-particle spectroscopic factors, where the latter cannot exceed the maximal natural 
occupation number fl24j, one can find that S ( a ] < >SL. 

A procedure for obtaining the TOF's on the basis of the TDM has been suggested in p3| . It is due to the particular 
asymptotic properties of the TOF's and is similar to the one suggested in pi| for deriving the one-body overlap 
functions from the one-body density matrix. 

In the case when two like nucleons (neutrons or protons) unbound to the rest of the system are simultaneously 
transferred, the following hyperspherical type of asymptotics is valid for the two-body overlap functions pi], [35], [56| 



*(r, R)-^N exp j (i? 2 + J^j (# + ir 2 )^ , (8) 

where r and R are the magnitudes of the relative and center-of-mass (CM) coordinates, r = ri r 2 and R = (ri +r 2 )/2, 
respectively, m is the nucleon mass and E = E^ — E^ is the two-nucleon separation energy. 
For a target nucleus with Jf ar — + the TOF in Eq. (||) can be written in the form 

® v jm{xx,x 2 ) = ^2{^> v jLs(ri,r 2 )'S>xs(cri,(T 2 )}j M , (9) 

LS 

where v is the number of the state of the residual nucleus with a given total momentum J, 

1 1 



XSM S (<TU 0-2) = {xiOl) ®Xi ( (72 )} SM = ( 2 ms i 2™ S! l 5Ms ) X 5™n ( ffl )^|™- 2 t^)' ( 10 ) 

m s1 m s 



and &vjlsm l (ri, 1*2) is the spatially dependent part of the overlap function. Performing a decomposition into angular 
momenta 1 = l r and L# (L — 1 + L^) corresponding to the relative and CM coordinates one obtains: 

$ v jSLMA*,K)=Y,$„ JSULR (r > R)\Y LR (R^®Y l (r)} . (11) 

fr* 1 L J LM L 



Then the TDM can be written as 

p^(x 1 ,x 2 ;x' 1 ,x' 2 ) = J2J212 PJSLIL R (r, R; r', R') A J s %l R (<T 1 ,a 2 ;r, R) o£; R'), (12) 

JM LS ILr S'L'l'L' R 
L'S' i>L' 

R 

where the radial part of the density matrix is 

P { ?SUL R (r, Rl r', R') = J2 K.isLiL R (r, R)$„js'm'L R (r' , R') (13) 

S' L' V L v 

and the spin-angular function is 

A J s ^ LR (a 1 ,a 2 ;r,R) = \ \y Lr (R) ®Y l (r)\ ® XsMsfa, ^)\ • (14) 
We will consider the diagonal part of the radial TDM in Eq. ( |l3| ) : 

PjsuL R (r,R-y,R') =J2®*»JSLiL R (r,m»JSLiL R (r',R'). (15) 
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For large r' = a and R' — b a single term with i/q, corresponding to the smallest two-nucleon separation energy, 



will dominate the sum in the right-hand side of Eq. (15). Then, according to Eq. (|8j), the radial part of the TOF 



&v JSLiL R (r, R) can be expressed in terms of the TDM as 

<2 Vo jSLiL R {r,H) = — 7 — r— -. t — -, (lb) 

®»oJSUL R {a,b) Nexp[-k^/(b 2 + \a 2 )j (b 2 + ±a 2 ) 7 

where k = (4m|i?|/?i 2 ) 1 / 2 is constrained by the experimental values of the two-nucleon separation energy E. 

The relationship obtained in Eq. (Ftq) makes it possible to extract TOF's with quantum numbers JSLILr from a 

— (21 

given TDM. The coefficient N and the constant k can be determined from the asymptotics of pjg LlL (r,R;r,R). 

III. RESULTS 
A. The two-proton overlap functions 

The procedure described in Section II has been applied to calculate the two-proton overlap functions in the 16 O 
nucleus for the transition to the + ground state of 14 C. The TDM obtained in |B3| in the framework of the low- 
order approximation (LOA) of the Jastrow correlation method has been used. In|33] the latter incorporates the 
nucleon-nucleon SRC in terms of the wave- function ansatz |p^7| : 

& A \r 1 ,r 2 ,...,r A ) = (C A )- 1 / 2 JT /(| r, - r,- \)^ D (n, r 2 , . . . , r A ), (17) 

l<i<j<A 

where Ca is a normalization constant and is a single Slater determinant wave function built from harmonic 

oscillator (HO) single-particle wave functions which depend on the oscillator parameter a osc ., having the same value 
for both protons and neutrons. Only central correlations are included in the correlation factor f(r), which is state- 
independent and has a simple Gaussian form 

f{r) = 1 - c exp(-/3V), (18) 

where the correlation parameter (3 determines the healing distance and the parameter c accounts for the strength of 
the SRC. The LOA keeps all terms up to the second order in h = f — 1 and the first order in g = f 2 — 1 in such a 
way that the normalization of the density matrices is ensured order by order p§[ ]. 

The values of parameters a osc . and j3 have been obtained |59| phenomenologically by fitting the experimental 
elastic formfactor data for 4 He, 16 and 40 Ca nuclei. The value of the parameter c has been determined |33| under 
the additional condition the relative pair density distribution 

p( 2 '(r)=|p (2) (r,R;r,R)dR (19) 

to reproduce at r = the associated value obtained within the Variational Monte-Carlo approach jlO|]. Thus, in 
the present calculations the following values of the parameters are used for 16 0: a osc . — 0.61 fm^ 1 , — 1.30 fm _1 , 
c = 0.77. 

In order to obtain the radial part of the TDM, p( 2 \r, R; r' , R') of Eq. (]l|), we use the analytical expression for 
the TDM obtained in substituting the coordinates of the two particles, ri and r2, by the CM R and relative r 
coordinates. Then, the TDM is multiplied by Agff L „ lllL ,, (eri, 02; r, R)Affil">i'"L'" (°i> °2> ^> R ') ana ^ the integration 
over the angles and summation over the spin variables lead to the radial part Pjs"L"i"L" ( r > ^! r '> ^')- 

S"'L"'l"'L'£ 

In order to obtain the values of the parameters k and N in Eq. ( |l6| ) simultaneously, we look for such a radial contri- 
bution P^jsLiL R ( r i R '> r '' R ') = K jSLiL R ( r > R )^u jSLiL R {r' , R') whose diagonal part p { ^ JSLlLR (r, R) minimizes the 
trace 



Tr 



P { JsLlL R ( r ' R ) ~ PvoJSLlL R ( r i R )) 



(20) 



The correct determination of these parameters requires a proper definition of the asymptotic region where the trace 
in Eq. ( |20j ) has to be minimized. If we denote the point in which P^j^slil r {t^ -^0 ^ maximum with (r m(1;c , 

Rmax), 
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FIG. 1: The 1 So two-proton overlap functions for the nucleus ls O leading to the + ground state of 14 C extracted from the 
JCM (left) and uncorrelated (right) two-body density matrices. 




FIG. 2: The 3 Pi two-proton overlap functions for the nucleus 16 leading to the + ground state of 14 C extracted from the 
JCM (left) and uncorrelated (right) two-body density matrices. 

the starting point of the asymptotic region (ro,i£o) is obtained looking for a point ro, at R = R m ax, for which 

(2) (2) 

Pjslil ( r o,Rmax) < 10% of p JSUL (r max ,R max ). When ro has been determined, we look for a point Ro, at 

(2) (2) 

r = ro, for which PjsLiL R ( r o> -^o) ^ 10% of pjg LlLR (ro,R max ). The length of the asymptotic region over r and R 
is determined by the requirement to obtain the separation energy which is maximally close to the experimental one. 
The asymptotic point (a, b) is chosen to be that one which gives the minimal least-squared deviation expressed by 
Eq. ©. 

When all the parameters are determined, Eq. ( fl6| ) can be used to calculate the radial part ^u jSLiL R (r, R)- Then, 
including also the spin-angular part in Eqs. (^TJ) and @ we obtain the TOF's. 

For a given set of quantum numbers JSLILr the TOF is calculated by minimizing the trace of the corresponding 
part of the TDM. Thus, for a particular final state J ir of the residual nucleus, different TOF's can be independently 
calculated using this procedure for each set of quantum numbers, and each one of them is fully responsible for the 
two-proton knockout process and the transition to the state J 7r . 

The TOF's obtained in the JCM for the 1 5o and 3 Pi states are presented in Figs, [j] and |[ respectively. They are 
compared with the uncorrelated TOF's obtained applying the same procedure to the uncorrelated TDM, i.e. with 
c = in the correlation factor of Eq. (|l8|). The notation for the partial waves in our case is 2S+1 Il- It differs from 
the generally accepted one 2S+1 lj because we have a different coupling scheme of spin and angular momenta. 

SRC depend on both relative and CM coordinates and it can be seen from the figures that they affect both size 
and shape of the TOF's. Their role, however, is different in the two states and, as it was already found in previous 
and different calculations (see, e.g., j?], ^2|), are much more important when the two protons are in the 1 Sq than in 
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the 3 Pi state. 

The spectroscopic factors corresponding to the 1 S and 3 Pi overlap functions are 0.958 and 0.957, respectively. 
Also the D wave can contribute for the transition to the + ground state of 14 C, but the corresponding TOF is very 
small and is not considered in the present study. 

As a next step, we derive the total TOF $„ji,|(i, X) in terms of a sum over different partial components, i.e. 

$ vJM (x,X)= J2 ^uJSUL R (r,R)A J s f lLR (a 1 ,a 2 ;r,R). (21) 

LSILr 

We integrate the squared modulus of the total TOF in Eq. (|2l|) over the angles and sum over the spin variables. 
The result can be written in the form (for the smallest value of v — vq): 



\$j M (x,X)\2=\$ JM (r,R)\ 2 = PjL8iL R (T,R), (22) 

LSIL 



where the bar denotes the integration over the angles and summation over the spin variables, and < &jm(^, R) is the 
radial part of the total TOF obtained after the integration and summation. Using the asymptotics of $ja/ {t, R) at 
r — > a, R — > b one can write: 



LSlLf 



3>jm{t,R) = f ; ^ tt 2 ■ (23) 

Nexp[-k^(b 2 + la 2 )}(b 2 + la2) 5/2 

The parameters N , k, a, b in Eq. ( |23| ) can be redetermined from the asymptotics of X)lszl.r P^jLSiL R ( r ' Q r i usm S 
the procedure already explained in the first part of this Section. Then, each partial radial component &jSLiL R (r, R) 
in Eq. (|2l|) can be separately calculated from Eq. ( |l6| ) using for each one of them the same coefficients N, k, a, b which 
correspond to the asymptotics of the total TOF. The asymptotic point (a, 6) determines the individual contribution 
of each partial overlap function to the total TOF. This prescription allows us to combine, with some approximations, 
the different radial components in Eq. (|2l|). 

The results for the 1 So and 3 Pi partial components have a similar behaviour as in Figs. [l| and g, the main difference 
is that they are somewhat reduced in magnitude. The reduction is determined by the contribution of each component 
to the total TOF. The spectroscopic factor corresponding to the total TOF is equal to unity in the uncorrelated case 
and 0.965 in the Jastrow case. 

B. The 16 0(e,e'pp) 14 C g . s . reaction 

The TOF's obtained from the TDM within the Jastrow correlation method have been used to calculate the cross 
section of the 16 0(e, e'pp) 14 C g . s . knockout reaction. 

Calculations have been performed within the theoretical framework of ||, q, 0. In this model the nuclear current 
operator is the sum of a one-body and a two-body part. The one-body part contains a Coulomb, a convective and a 
spin term. For pp knockout the two-body current contains only the contributions of non charge- exchange processes 
with intermediate A-isobar configurations (t], [l4|] . In the final state the mutual interaction between the two outgoing 
nucleons is neglected and the scattering state is given by the product of two uncoupled single-particle distorted wave 
functions, eigenfunctions of a complex phenomenological optical potential which contains a central, a Coulomb and a 
spin-orbit term Jjlj . 

Numerical results are shown in Fi gs. [3 and ^ for two kinematical settings considered in the experiments performed 
at NIKHEF pi [l7|, [H| and MAMI pa" 2^. In Fig. g the cross section is calculated in the super-parallel kinematics 
of the MAMI experiment, where the two nucleons are ejected parallel and anti-parallel to the momentum transfer 
and, for a fixed value of the energy to and momentum transfer q, it is possible to explore, for different values of the 
kinetic energies of the outgoing nucleons, all the possible values of the recoil momentum pe of the residual nucleus. 
In the calculations the incident electron energy is fixed at Eq = 855 MeV, uj = 215 MeV and q = 316 MeV/c. In 
Fig. ^ a specific kinematical setting included in the experiments carried out at NIKHEF is considered, with Eq = 584 
MeV, w = 212 MeV and q = 300 MeV/c. The kinetic energy of the first outgoing proton T{ is 137 MeV and the 
angle 71, between the outgoing proton and q, is 30° on the opposite side of the outgoing electron with respect to 
the momentum transfer. Changing the angle 72 on the other side, different values of the recoil momentum pb are 
explored in the range between —250 and 300 MeV/c, including the zero values at 72 — 120°. 
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FIG. 3: The differential cross section of the 16 0(e, e'pp) 14 C g . s . reaction as a function of the recoil momentum ps in the 
superparallel kinematics with Eo = 855 MeV, uj = 215 MeV and q — 316 MeV/c. Positive (negative) values of pb refer to 
situations where pb is parallel (anti-parallel) to q. The curves are obtained with different treatments of the TOF: 1 So (dashed 
line) and 3 Pi (dotted line) as independent TOF's in the JCM in the left panel and as partial components in the right panel, 
the total TOF (solid line), the TOF from the spectral function (SF) jj], (dot-dashed line), the product of a pair function of 
the shell model and the correlation function of Eq. (jig) (SM+CORR) (dot-dot-dashed line). 
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FIG. 4: The differential cross section of the 16 0(e, e'pp) 14 C g . s . reaction as a function of the angle 72 in a NIKHEF kinematics 
with E = 584 MeV, w = 212 MeV, q = 300 MeV/c, T{ = 137 MeV and _7i = -30°, on the opposite side of the outgoing 
electron with respect to the momentum transfer. Line convention as in Fig. |3[ 



The cross sections calculated with the 1 S$ and 3 Pi TOF's as independent and fully responsible for the knockout 
process are displayed in the left panels of Figs. || and |[ In the right panels the cross sections obtained with the 
total TOF from the Jastrow TDM are plotted and compared with the contributions given by the 1 S'o and 3 Pi partial 
components. 

These results are compared in the figures with the cross sections already shown in j?]], where the TOF is taken 
from a calculation of the two-proton spectral function (SF) where a two-step procedure has been adopted to 
include both SRC and LRC. LRC are calculated in a shell- model space large enough to incorporate the corresponding 
collective features which influence the pair removal amplitude. The single-particle propagators used for this dressed 
random phase approximation (RPA) description of the two-particle propagator also include the effect of both LRC 
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and SRC. In the second step that part of the pair removal amplitudes, which describes the relative motion of the pair, 
is supplemented by defect functions obtained from the same G-matrix which is also used as the effective interaction in 
the RPA calculation. Different defect functions are produced by different realistic NN potentials. The results shown 
in Figs. H and || are obtained with the Bonn-A potential. The explicit expression of the TOF's is given in a form of 
the same kind as in Eq. (|2l|), in terms of a combination of CM and relative wave functions. The 1 5'o and 3 Pi relative 
waves give the main contribution for the transition to the + ground state, while only a negligible contribution is 
given by the D wave. The results of this model are able to give a proper description of available data |l7], [H| . 

In the figures are also shown for a comparison the results obtained with a simpler approach, where the two-nuclcon 
wave function is given by the product of the pair function of the shell model and of a Jastrow type central and 
state independent correlation function. In this approach (SM+CORR) the ground state of 14 C is described as a pure 
(lpi)~ 2 hole in 16 0. In order to allow a more direct and clear comparison with the TOF's from the Jastrow TDM, 
HO single-particle wave functions and the same correlation function as in Eq. (|l^), with the same parameters as in 
the calculation of the TDM, have been adopted. 

The shape of the calculated cross sections is determined by the value of the CM orbital angular momentum Lr, 
that is — for 1 Sq and Lr = 1 for 3 Pi. When the two components are combined in the TOF the shape is driven 
by the component which gives the major contribution, that is Lr — and So at lower values and Lr = 1 and 3 Pi 
at higher values of the recoil momentum. 

The role of correlations and two-body currents is different in different relative states. SRC are quite strong and even 
dominant for the 1 S'o state and much weaker for the 3 Pi state. Moreover, the role of the isobar current is strongly 
suppressed for 1 So pp knockout, since there the generally dominant contribution of that current, due to the magnetic 
dipole NN <-> NA transition, is suppressed |H| ^2|. Thus, the role of SRC is emphasized in 1 Sq knockout, while the 
A current is emphasized in 3 Pi knockout. These general features, which have been found in all the previous studies 
of the exclusive 16 0(e, e'pp) 14 C reactions, are confirmed, for both kinematical settings here considered, also in the 
present calculations. The ^-Sq results shown in the figures are dominated by the one-body current and thus by SRC, 
while the A current gives the main contribution to the 3 Pi results. 

One of the main results of the previous theoretical investigations is the dominance of 1 Sq pp knockout in the 
16 0(e, e'pp) 1 Cg.s. reaction. These theoretical predictions have been clearly confirmed in comparison with data. Even 
though the contribution of 3 Pi pp knockout can become important and even dominant at large values of pb, it is clear 
that a TOF where only the 3 Pi state is included is unable to give a reliable description of the two-proton knockout 
process. In contrast, it can be seen from Fig. |^ that the cross section calculated with the Jastrow TOF for the 1 S'o 
state is close to the SF and also to the SM+CORR results at low values of pb, up to ~ 150 — 200 MeV/c, that is 
just in the region where the 1 S'o contribution is dominant. For pb > 200 MeV/c 3 Pi knockout becomes dominant 
with all the different treatments of the TOF. The results with the 3 Pi TOF from the Jastrow TDM is however much 
larger than the SF result and also larger than the SM+CORR cross section. It can be noted that even the 1 5o curve 
in Fig. H is, at large values of the momentum, higher that the SF result. This is an indication that SRC in the JCM 
produce a stronger enhancement of the high-momentum components. 

The behaviour of the pure 1 S'o result in the left panel of Fig. [| is somewhat similar to that of the SF and SM+CORR 
cross sections, which appear driven by the contribution. There are anyhow significant differences in the shape 
and large differences in the size of the various results. 

The cross sections calculated with the total TOF, where the 1 5'o and 3 Pi partial components from the Jastrow 
TDM are combined, are shown in the right panels of Figs. || and |[ It can be seen that in both kinematical settings 
the 1 Sq component dominates at low values of pb, while the 3 Pi component produces a strong enhancement of the 
cross sections at high momenta. The contribution of the partial ^^Sq component is reduced with respect to the results 
in the left panels, where 1 S'o is fully responsible for the knockout process. Thus, the cross sections calculated with the 
total TOF from the Jastrow TDM are somewhat reduced at low recoil momenta. Also the contribution of the partial 
3 Pi component is slightly reduced with the respect to the 3 Pi results displayed in the left panels. The contribution 
of the 3 Pi component to the total TOF is, however, much more relevant than with the other theoretical treatments 
considered in the figures and the enhancement at high momenta turns out to be much larger. Thus, the shape of the 
cross sections with the total TOF from the JCM is flatter than with the SF and SM+CORR results. 

In comparison with the SF calculations, the cross sections with the Jastrow TOF are lower at low recoil momenta 
and much larger at high momenta, due to the larger contribution of the 3 Pi component in the TOF. The SM+CORR 
cross sections are higher than the other results at low recoil momenta. This is an indication of a stronger contribution of 
SRC in this calculation. This contribution, however, depends on the particular expression adopted for the correlation 
functions, that in the calculations of Figs. || and [| is exactly the same as in the calculation of the TDM. At high 
momenta the SM+CORR cross sections remain always higher than the SF results, but generally lower than the results 
given by the TOF from the JCM. 

Although obtained from a calculation of the TDM within the JCM where only SRC are included the TOF used 
in our calculations are able to reproduce the main qualitative features of the 16 0(e, e'pp) 14 C g . s . cross sections which 
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were found in previous theoretical investigations and also in the analysis of the available data. This means that 
the procedure suggested in Q to calculate the TOF's from the TDM can be applied and exploited in the study of 
two-nucleon knockout reactions. 

The large differences found in Figs. || and || indicate that the calculated cross sections are very sensitive to the 
different approaches used and to the theoretical treatment of nuclear structure and correlations in the TOF. It would 
be interesting to apply the procedure used in this work for the calculation of the TOF's to more refined treatments 
of the TDM. 



IV. CONCLUSIONS 

The results of the present work can be summarized as follows: 

i) The two-nucleon overlap functions (and their norms, the spectroscopic factors) corresponding to the knockout of 
two protons from the ground state of 16 O and the transition to the ground state of 14 C are calculated using the 
recently established relationship between the TOF's and the TDM. In the calculations the TDM obtained 
within the JCM |53) is used. Though only SRC are accounted for in the Jastrow TDM, the results can be 
considered as a first attempt to use an approach which fulfils the general necessity the TOF's to be extracted 
from theoretically calculated TDM's corresponding to realistic wave functions of the nuclear states. Of course, 
the quality of the results will depend heavily on the availability of a realistic TDM incorporating all necessary 
types of NN correlations. 

ii) The contributions of the two-proton overlap functions corresponding to the removal of 1 Sq and 3 Pi pp-pairs 
from 16 are calculated in two manners: 1) when each one is fully responsible for the knockout process, and 
2) when they are partial components of the total TOF. The 1 5'o and 3 Pi results obtained in the two manners 
are similar, the main difference being that the partial components in case 2) are reduced in magnitude. The 
comparison between the results of correlated (Jastrow) and uncorrelated TOF's shows that SRC depend on 
both relative and CM coordinates and affect the size and the shape of the 1 5'o and 3 Pi overlap functions. The 
effects of SRC, however, are much stronger when the two protons are in an 1 Sq state. 

iii) The TOF's obtained from the Jastrow TDM are included in the theoretical approach of |3|, ^| Q to calculate the 
cross section of the 16 0(e, e'pp) 14 C g . s . knockout reaction. Calculations are performed in two different kinematics 
that have been realized for the cross section measurements at NIKHEF and MAMI. The results are compared 
with the cross section calculated, within the same theoretical model for the reaction mechanism, with different 
treatments of the TOF, in particular with the TOF obtained from a calculation of the two-proton spectral 
function of 16 O [in |2^] where both SRC and LRC are included. The calculated cross sections are very sensitive to 
the theoretical treatment and different results are produced by the different TOF's. The cross sections calculated 
in the present work, where the TOF's are extracted from the Jastrow TDM, confirm the dominant contribution 
of 1 S'o pp knockout at low values of recoil momentum up to ~ 150 — 200 MeV/c. The 3 Pi contribution is mainly 
responsible for the high-momentum part of the cross section at pb > 200 MeV/c. 

iv) Our method is applied in the present work only to the ground state of 14 C. It can be used also for the excited 
states. Our main aim was to check the practical application of all steps of the method for a given state of 
the residual nucleus. Therefore, the results obtained for the 16 0(e, e'pp) 14 C g . s . reaction, which are able to 
reproduce the main qualitative features of the experimental data and of the cross sections calculated with 
different treatments of the TOF's, can serve as an indication of the reliability of the method, that can be 
applied in a wider range of situations and to more refined approaches of the TDM. 
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